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ABSTRACT 



Over the last few years, needlets have emerged as a useful tool for the analysis 
of Cosmic Microwave Background (CMB) data. Our aim in this paper is first 
to introduce in the CMB literature a different form of needlets, known as Mexi- 
can needlets, first discussed in the mathematical literature by Geller and Mayeli 
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(2009a,b). We then proceed with an extensive study of the properties of both 
standard and Mexican needlets; these properties depend on some parameters 
which can be tuned in order to optimize the performance for a given apphcation. 
Our second aim in this paper is then to give practical advice on how to adjust 
these parameters for WMAP and Planck data in order to achieve the best prop- 
erties for a given problem in CMB data analysis. In particular we investigate 
localization properties in real and harmonic space and propose a recipe on how 
to quantify the influence of galactic and point source masks on the needlet co- 
efficients. We also show that for certain parameter values, the Mexican needlets 
provide a close approximation to the Spherical Mexican Hat Wavelets (whence 
their name), with some advantages concerning their numerical implementation 
and the derivation of their statistical properties. 

Subject headings: (cosmology:) cosmic microwave background — cosmology: 
observations — methods: data analysis — methods: statistical 



Introduction 



Over the last decade, wavelet systems have grown as one of the most important tools 
in the analysis of Cosmological and Astrophysical data. A lot of proposals for wavelet 



systems on the spher e have been advanced in the mathematic a. 



mm 



Rosea 



mm . 



Freeden fc Schreiner I ( 1998).lAntoine fc Vanderghevnst ( 2007 ). 



Wiaux et al. I (120071 ) . IStarck et al. I (120061 ) and the references therein. 



literature, see for instance 



Dahlke et al.l ( 120071 ). iHolschneider fc IdeA^ 



Some of these a ttempts have be e n exp licitly motivated by Astronomy and/or Cosmology 
(see for instance iMcEwen et al. I (120071 ) for a review). In particular in the area of Cosmic 
Microwave Background (CMB) data analysis, wavelets have been used for a large number 
of applications (see references in the next paragraph). The interest for wavelets in this area 
is very easily understood; predictions from CMB theory are typically cast in the Fourier 
domain, however exact Fourier analysis cannot be entertained because of the presence of 
foreground and masked regions. The double-localization properties of wavelet systems (in 
real and harmonic domain) hence turn out to be most valuable. Moreover, addressing im- 
portant issues such as the possible existence of features and asymmetries in CMB maps is 
nearly unfeasible without ideas which are broadly related to the wavelet literature. 

Among spherical wavelets, particular attention has been recently de voted to so-called 



needlet s, which were introduced into the Functional Analysi s literature bv iNarcowich et al. 



( l2006al Jbl): their statistical properties were first considered by lBaldi et al. I ( l2009al Jbl). Needlets 
enjoy several properties that make them worth of attention for Cosmological data analysis. 
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In particular, they are cor nputationally very simple, and naturally adapted to standard 



packages such as HealPix ( IGorski et al.l l2005l ): they do not require any form of tangent 



plane approximation, but they are naturally embedded into the manifold structure of the 
sphere; they are compactly supported in the harmonic domain, i.e. they depend only on a 
finite number of multipoles which are explicitly known and can be controlled by the data 
analysts; they are quasi-exponentially localized in real space, i.e. their tail s decay faster 



than any polynomial; and finally, it has been shown in iBaldi et al. I ( l2009al ) that random 
needlet coefficients enjoy a very useful uncorrelation property: namely, for any fixed angu- 
lar distance, random needlets coefficients are asymptotically uncorrelated as the frequency 
parameter grows larger and larger. As well-known, uncorrelation entails independence in 
the Gaussian case: as a consequence, from the above-mentioned property it follows that 
needlet coefficients from a CMB map can be seen as nearly independent at high frequencies, 
making thus possible the introduction of a variety of statistical procedures for testing non- 
Gaussianity, estimating the angular power spectrum, testing for asymmetries, implement- 



ture data, and many others, see for instance 


Baldi et al. 


f 2009a 


b). 


Gi 


Liilloux et al. 


(2009 


), 


Lan & Marinucci 


f 2008b 


), 


Pietrobon et al. 


(2006), 


Marinucci et al. 


(2008). 


Fav, G. et al. 


(2008 




Rudjord et al. 


( 


2009a 


lb), 


Cabella et al. 


OT). 


More recently, the needlet con- 



struction has also been e xtended to the case of s pin /p olarization data, see for i nstance 
Geller et al. I J2008[ l2009l ). iGeller fc Maveh I J2009dl ) and lGeller fc Marinucci I f l2008h . 



The first purpose of this paper is to introduce a new kind of needlets to the field of CMB 
analysis following an appr oach which has been very recently advocated in mathematics by 
Geller fc Mayeli I ( 2009aU blJd). This approach (which we shall discuss in Section 2) can be 



labeled Mexican needlets. As we shall discuss below, a special case of the Mexican needlets 
provides at high frequencies a close app roximation to the ^y i dely used Spherical M exican 



Hat Wavelet (SMHW, see for instance ICavon et al. I (120011 ): IVielva et al. I (120041 ) with 



some advantages in terms of their numerical implementation and the investigation of their 
localization and statistical properties. As such, the investigation of their properties in this 
case will allow to understand the stochastic properties of SMHW and compare them with 
standard needlets. Mexican needlets depend on a parameter p, and we shall show how this 
parameter can be tuned to improve the localization properties in the real or in the harmonic 
domain. 

The second purpose of this paper is to provide a practical description of properties of 
different needlet types important for CMB analysis. Proper knowledge of the localization 
properties on the pixelized sphere as well as in multipole space is crucial for selecting and 
applying the proper type of needlet to a specific problem. Although the exact mathematical 
properties of the needlets are well known, we have studied the properties which are of high 
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importance for CMB analysis and which are too comphcated to be easily deduced from the 
mathematical results. In particular, in the presence of foregrounds and masks, it is important 
to know their influence on the needlet coefficients. For the SMHW it has been shown for 



several application s (IVielva et al. I |2004J : iMukherjee &: Wang I l2004l : iMcEwen et al. 1 12004 



Curto et al. I l2009l ) that an extended scale dependent mask must be used when analyzing 
masked data with wavelets. Here we will study this in detail for the different needlet types 
with different parameter values. 

We then provide a very thorough comparison between different needlets. The previous 
discussion leads naturally to the issue about their optimal construction, i.e. how to devise 
numerical recipes which will enhance their localization p r operties. Here, we shall compare 



the numerical recipe implemented by lBaldi et al. I (l2009al ): iPietrobon et al. I (|2006[) w ith an 



alternative proposal based on Bernstein polynomials (see also iGuilloux et al. I ( 120091 ) for a 



related numerical investigation). The latter entail weight functions with just a finite number 
of bounded derivatives; to distinguish it from the previous construction we will label this 
procedure Bernstein needlets. We stress, however, tha t the underlying mathem atical theory 
presents no real novelty as compared to the results by iNarcowich et al. I (j2006aU bl) . We then 
go on to provided vast numerical evidence on the various forms of localization, by means of 
a number of different indicators. In particular, the role of the different parameters in the 
determination of the various properties is fully exploited. 

As well-known, there is usually a trade-off between localization properties in the fre- 
quency and real domains, as a consequence of the Uncertainty Principle ( "It is impossible 
for a no n-zero function and i t s Fou r ier transform to b e simu ltaneously very small", see for 
instance INarcowich fc Ward I (119961 ): iHavin fc Joricke I (119941 )): the main purpose of this pa- 



per is to show how Mexican and standard needlets jointly provide a flexible set of tools 
where each user can optimize this trade-off according to the needs of a speciflc data analysis 
problem. 

The plan of the paper is as follows: in Section |2l we review the standard needlet con- 
struction, introduce Mexican needlets and compare their respective properties from a math- 
ematical point of view. In the following sections, we provide numerical evidence on the 
localization properties of these procedures by means of several different flgures of merit; 
we discuss at length the interplay among the different properties and the trade-off to face 
when choosing which procedure to adopt for a given Astrophysical problem. In section [TJ we 
summarize the main indicators used and the properties of the different needlets measured in 
terms of these indicators. In the Appendix, we provide some details on the numerical recipes 
we adopted, some short discussion on mathematical properties and analytic fits. 
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Mexican and Standard Needlets 



Th e construction of the stan dard needlet system is detailed in lNarcowich et al. 



see also Marinucci et al. 



J2006al lbh. 



( 120081 ): we sketch here a few details to fix notation and we provide 
in the Appendix a more detailed dis c ussion for completen ess. The introd u ction of Mexican 
needlets is due to Geller fc Mayeli ( 2009a b^) (see also Freeden et al. ( 1998 )): they are 
used here for the first time in the Astrophysical literature, so we provide below a more com- 
plete discussion. We point out that for p = 1, an analogous proposal was first advocated in 



the Astrophysical literature by ISanz et al. I ( l2006l ). 



The basic needlet function can be described in real space as follows: 



em(x) 



Here, x refers to a position {6, 0) on the sphere, Yim are spherical harmonic functions, j is 
the scale (frequency) of the needlet and {Xjk} is a set of cubature weights corresponding to 
the cubature points {C,jk} ] for simplicit y, they can be tak en to be equal to the pixel areas 
and the pixel centres in the HEALPix ( iGorski et al.l 120051 ) grid used for CMB analysis, i.e. 
we shall consider Xj^ = Xj = An/Nj, where Nj is the number of pixel in the pixelization 
we are working with. The needlet function itself is contained in the function bi{B,j) (or 
bi for short) in harmonic space, B being one of the parameters deciding the properties of 
the needlet. The difference between the needlet systems we are going to discuss can thus be 
traced in the form of the weight function 6^. 

1) Standard needlets: Let 0(^) be an infinitely differentiable (i.e., C°°) function 
supported in |^| < 1, such that < 0(0 < 1 and 0(0 = 1 if I0 <l/B, B >1. Define 



B' 



0(0 > SO that \/i>B , J2 ^^(;^) = 1 • 

j=0 



(2) 



For a 



For standard needlets we then obtain 6^ from this function 6(0 by bi{B,j) = b(-^ 
given scale j, the needlet function in harmonic space is centered at a multipole i* ~ B^ . Thus 
a given scale j is mainly influenced by multipoles close to i*. It is immediate to verify that 
6(0 7^ only if < |0 ^ B. An explicit recipe t o construct a funct i on 6(0 with the previous 



features is discussed in Appendix |A] (compare iPietrobon et al. I (120061 ) . iMarinucci et al. 
(120081)') . The main localization property of needlets is established in iNarcowich et al. 



(l2006aj), where it is shown that for any M G N there exists a constant cm > s.t., for every 
e G S^: 
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where d{^jk,^) denotes the usual distance on the sphere. More exphcitly, needlets are almost 
exponentially localized around any cubature point, which motivates their name. 

2) Bernstein Needlets The bound which we just provided to establish the localization 
properties of needlets depends on some constants cm which we did not write down explicitly. 
Such constants depend on the form of the function b{^), and turn out to be rather large in the 
case of standard needlets. In Appendix [B], we give another method of construction of b{^), 
where such function is no longer infinitely differentiable but rather has a finite num ber of 



bounded derivatives. The localization theory described in iNarcowich et al. I fcoOGallbh 



goes 



through without any modification, as do the stochastic properties established by lBaldi et al. 
(l2009a| ). We do no longer have quasi-exponentially decaying tails, however, but it is pos- 
sible to establish a weaker result, namely the decay with a polynomial rate, depending on 
the number of bounded derivatives we are allowing for b{C,). It may hence seem that this 
construction should enjoy worse properties - but in practice this is not the shown by 

our simulations in the sections to follow. As for standard needlets, we have be{B,j) = b{^), 
but note that is different for Bernstein needlets as detailed in the Appendix. As for the 
standard needlets, the needlet function in harmonic space is centered at £* ~ . 



3) Mexican Needlets The construction in lGeller fc Mayeli I (l2009a( ) is similar to stan- 
dard needlets, insofar as a combination of Legendre polynomials with a smooth function is 
proposed; the main difference is that for standard needlets the kernel is taken to be com- 
pactly supported (i.e., depending only a finite number of multipoles i), while the Mexican 
needlet construction draws information from all frequencies at any scale. More precisely, we 
shall consider weight functions be{B,j) of the form 



HB,j) = {^r'e~^ , (3) 
for p = 1, 2, 3, ... For instance, for p = 1 the Mexican needlet takes the form 



and for higher p we have 

B'^i' 4:7T 



p2 Of _i_ ] 



Indeed, for mathematical rigour £^ should be replaced by the eigenvalue i{i + 1), but for 
CMB data analysis the difference is negligible and we shall use i"^ for notational simplicity. 
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As mentioned before, Mexican needlets are not supported on a finite number of mul- 
tipoles, so the discussion of their locahzation properties in the harmonic domain requires 
some care. Moreover, because we need to focus on an infinite number of spherical harmon- 
ics, from a strictly mathematical point of view exact cubature a nd reconstruction formul ae 



cannot hold. Nevertheless, it must be added that the approach by lGeller fc Mayeli I (l2009al Jbl) 



enjoys some undeniable strong points, some of which we list as follows: 

1) Mexican needlets enjoy extremely good localization properties in the real domain; 
more precisely, at a fixed angular distance x their tails decay as exp{—B'^^x'^/4), as j grows 
to infinity. 

2) By adjusting the parameter p, one has available a family of wavelets which can be 
optimized in terms of the desired localization properties (as we shall show below, a growing 
p improves the localization in the harmonic domain and decreases the localization in the real 
domain) 

3) The previously mentioned mathematical issues on the cubature points are largely 
negligible from a numerical point of view 

4) The Monte Carlo evidence provided below proves that Mexican needlets compare 
favorably with standard needlets under a variety of circumstances and for many different 
indicators 

5) Analytic expressions can be provided for their high-frequency behavior in real space. 

Concerning the last point, it is important to remark the following. It can be shown 
that Mexican needlets for p = 1 provide a very close approximation of the widely popular 
Spherical Mexican Hat Wavelets (SMHW), see Appendix O Even in this case, though, the 
implementation through needlet ideas in our view yields important benefits: 

a) the weight function is explicitly given, making easier the implementation and the 
validation of numerical codes 

b) the localization structure in harmonic domain can be analytically studied and con- 
trolled 

c) the correlation structure of random Mexican needlet coefficients is explicitly given 
and can be used for statistical inference 

d) the range of scales to be considered to retain the information from the data is mathe- 
matically determined in terms of the frequencies j, rather than by an ad hoc choice of scales 
in the real domain ClS db function of angular distance. 
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3. Correlation properties of standard and Mexican needlets 

In the sequel of the paper, we shall compare three properties of these needlet construc- 
tions, namely their localization in the real domain, the localization in the harmonic domain, 
and the statistical properties of needlet coefficients, primarily their correlation structure. 
Spherical needlet coefficients are defined as 



'jk 



"^^^ " m=-l 



where T{x) is the CMB temperature field. The correlation coefficient is hence given by 

{(3jk(!^jk') 



Corr {I3jk,l3jk' 



_ Y.,>MB.3Y-^c,p,{d{i,k.i,k')) 

where Pe is the Legendre polynomial of degree i and Ce is the power sp ectrum of the CMB 



tempe rature field; the last step follows from the well-known identity ( iVarshalovich et al. 



(119881)) 



'L 2/ -I- 1 



Air 



For standard needlet coefficients, it was shown by iBaldi et al. I (|2009aJ) that under general 
conditions the following inequality holds 



\Corr {(3jk, (3jk')\ < 



M 



'1 + B^d{^,k,^jk')) 



M ' 



some C^j > , 



(6) 



where d{^jk,C,jk') is the standard geodesic distance on the sphere. In words, for any two 
points at a finite distance on the sphere the correlation between needlet coefficients centred 
on this points decays to zero as the frequencies grow larger and larger. Of course, under 
Gaussianity this simply implies that the coefficients become nearly independent at high 
frequencies. For Mexican needl e t coeffi c ients, th e situation is slightly more complicated, as 
discussed by iLan fc Marinucci I (l2008al ) ; iMayeli I (120081 ) . More precisely, let us assume that 
the CMB angular power spectrum behaves as 



Ce — {\aemf) 



Gi£)t 



(7) 
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where G{i) is some smooth function, for instance the ratio of two positive polynomials. 
Clearly ([7]) provides a good approximation to CMB spectra, with spectral index a ^ 2. It is 
then possible to show that, for p such that a < 4p + 2, there exist some constant Cm > 
such that 



Cm 

\Corr {I3jk;p, f3jk'-p)\ < 



(4p+2-a) 



where some po ssible logarithmic factors have been neglected, see iLan fc Marinucci I (j2008al ): 



Maveli I (|2008|) for details 



We should note that while uncorrelation holds for standard needlets no matter what 
the rate of decay of the angular power spectrum, here we need a not to be "too large" 
as compared to the order of the Mexican needlet we are using. The intuition behind this 
result is the following. Mexican needlets (and similarly Spherical Mexican Hat Wavelets) are 
not compactly supported in the harmonic domain; in other words, whatever the frequency 
j, they are drawing information from the smallest multipoles, i.e. those most affected by 
cosmic variance. The faster the decay of the spectrum (i.e., the higher the a), the greater the 
influence of this low frequency components on the behavior at high j. In order to compensate 
for this cosmic variance effect, it is necessary to ensure that the Mexican needlet filter will 
go to zero fast enough in the harmonic domain. Clearly, the higher the p, the faster our 
wavelet will approach zero at low multipoles, thus compensating for Cosmic Variance effects. 
In practice, however, for CMB data as we mentioned before a can be taken to be equal to 
2, whence the correlation coefficient is seen to decay to zero even for the smallest p = 1. 
Indeed, our numerical results below will show that for physically realistic angular power 
spectra Mexican needlets outperform standard ones in terms of uncorrelation properties, 
thus providing one more possible motivation for their use on CMB data. 

To close the introduction to Mexican and standard needlets we present plots showing 
a comparison between the Mexican needlets for p = 1 and SMHW (see figures [T]), Mexican 
needlets for different values of the parameter p at different frequencies j (see figures [2]) and 
the weight function be = b£{B,j) for Mexican and standard needlets (see the 2 top figures 
in figures [3] for comparison between Mexican and standard needlets, see bottom figure in 
figures [3] for seeing how the Mexican weight function depend on the parameter p). 

Note that for spherical Mexican needlets, the needlet function in harmonic space is 
no more centered at i* ~ since one no longer has a symmetric distribution around the 
maximum of 6^, as can be seen in figures [3] . To obtain a measure of the multipole we are 
looking at, we introduce therefore a weighted average defined the following way: 
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t U, B) 



which will be used for all kinds of needlets. 

Comparison: Mexican needlets and SIVIHW at low multipoles 




SMHW,R= 35.19[deg] 
mex,p=1,j=1 
SMHW,R=23.16[deg] 
mex,p=1 ,j=2 
SMHW,R=14.59[deg] 



mex,p=1 ,j=3 



150 



e[deg] 



Comparison: Mexican needlets and SMHW at higli multipoles 



500 



400 



300 



> 200 



100 



-100 
0, 



SMHW,R= 0.340[deg] 
mex,p=1 ,i=1 1 
SMHW,R=0.212[deg] 
mex,p=1J=12 
SMHW,R=0.0830[deg] 
mex,p=1 ,j=14 



0.2 



0.4 



0.6 
e[cleg] 



0.8 



1.0 



1.2 



(9) 



Fig. 1. — Comparison between Mexican needlets (grey) and SMHW (black). On the top we 
present the results at low multipoles I (big angular scales) , while on the bottom the results 
for high multipoles I (small angular scales) . As expected the higher the multipoles the better 
the similarity. 
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Mexican needlet j=1 Mexican neediet j=6 




0.0 0.5 1.0 1.5 2.0 0.0 0.2 0.4 0.6 

e[deg) e[deg) 



Fig. 2. — Mexican Needlets for B = 1.6 for different values of p and j. Upper left plot: j=l, 
upper right plot: j=6, lower left plot: j=ll, lower right plot: j=16 

4. Real space localization: distance of influence from a mask 

For the apphcation of the needlet transform to CMB analysis, it is crucial to know 
how well the needlet coefficients are localized on the sphere. This is particularly true in the 
presence of foreground contaminants (diffuse galactic foregrounds or cxtragalactic sources) 
or in cases were parts of the CMB maps have been masked. For all practical analysis one 
will always need to ask the questions 'How far away from the galactic plane are the needlet 
coefficients unaffected by the galaxy?', or 'How far away from the mask are the needlet 
coefficients expected to behave as if the mask were not present?'. 

As there is an infinite number of different needlet bases, one cannot run simulations 
each time one changes the basis in order to infer the localization properties for this particular 
kind of needlet. By comparing the influence of the mask in simulations of some needlets to 
expected properties based of these needlet functions, it is our aim to obtain an understanding 
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Weight function for H/lexican and Standard at low frequencies j 



b„ mex, p=1 and j=1 
b„ std, j=1 

b„ mex, p=1 and j=4 
b„ std, j=4 
b„ mex, p=1 and j=6 
b„ std, j=6 




20 40 60 80 

I 

Weight function for mexican and standard at high frequencies j 




Weight function b, for mexican needlets for different vaiues of p 




b„ mex, p=1 and j=14 
b„ mex, p=2 and i=14- 
b„ mex, p=3andj=14- 
b„ mex, p=4andj=14- 



Fig. 3. — Mexican and standard needlet weight functions for B = 1.6 for different values 
of p and j. The upper two plots show the comparison between Mexican {p = 1) and standard 
needlets for low and high frequencies whereas the lower plot compares Mexican needlets for 
different values of p. 
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of the effect of the shape of the needlet functions to simulated CMB maps and thereby infer 
formulae which can be applied to a large group of needlets without the need of running new 
simulations each time. 

Our main goal is to find a relation for the minimum distance from a mask where the 
needlet coefficients are not significantly affected by it. We perform simulations and compare 
the needlet coefficients with and without the presence of the mask in order to define the size 
of the contaminated regions. We found that the starting point for the most stable way to 
define the contaminated and safe regions was to construct the correlation coefficient between 
the masked and unmasked needlet coefficients. The higher the correlation, the less the 
infiuence from the mask. We have obtained these correlation coefficients from an ensemble 
of 10000 simulations. Each simulations was created and treated according to the following 
procedure: 



1. Using the WMAP (IHinshaw et al.l l2009l ) best fit power spectrum C^, we generate a 
random set of harmonic coefficients aim and then transform to obtain the corresponding 
temperature map T{6,(f)). 

2. Make a needlet transform from the map T{6, 0) and obtain the needlet coefficients (at 
pixel k) where the superscript nm stands for "no mask"; 

3. Save the quantity in order to be able to be able to calculate the variance of the 
needlet coefficients at the end of the iterations; 

4. Multiply the temperature map T{6, 0) above with the mask and again make the needlet 
transform in order to obtain the coefficients (3Jl, where the superscript m stands for 
"with mask" ; 

5. Again, save the quantity {PJlY', 

6. Save also the square of the difference between the two maps, [P^^ — PjlY', 



After running all simulations, we are interested in constructing the correlation coefficient 



were () represents average over simulations and over all pixels k at the same distance 9 from 
the mask. CJ^'^"^{6) is then the correlation coefficient as a function only of the distance 6 
from the mask for a given scale j. 
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The fraction of influence from the mask at a given 9 is then given by 1 — Cj^'^"^{6). 
We will now define the critical angle dcrit which is the distance from the border of the mask 
after which the total fraction of the influence of the mask is smaller than a certain threshold 
r. Thus for a threshold of r = 0.01, 99% of the influence of the mask is found inside the 
critical angle. This area should then be masked if the remaining 1% of the total influence is 
accepted for the data to be used. The thresholds we choose for both masks lie in the range 
[0.1, 10~^]. In mathematical terms, the critical angle is defined as 

where is the border of the mask. Note that Ocrit is defined to be zero at the border of the 
mask so that it directly measures the critical distance from the border of the mask. 

In the following two subsections we will study the result of this procedure to obtain the 
critical angle applied to galactic and point source masks based on the WMAP and Planck 
experiments. 

Since these runs were computationally heavy, we chose to do them with Nside = 512 and 
^max = 1200 and only for the Planck point source holes we needed to increase to Ngi^e = 1024 
and ^max = 2100. This has the consequence that when making the needlet transform, where 
the respective harmonic decomposition coefficients aim are multiplied by hi for high values 
of j, it is an "incomplete transform", because at £ = 1200 respectively ^ = 2100, bi is not 
yet tending to zero for high values of j. The consequence is that for high £*'s the values of 
the critical angle start to grow. We discard the values of the critical angle for those £*'s. 
Additionally, to have a better estimate of i* we calculated it (via equation [9]) using the 
functions bi for imax = 45OC0. 

The calculations of the correlations C™'"™'(^) were done for: 4 (point source holes) and 
5 (galactic cut) values of the threshold r, three values of i? = 1.6, 2 and 1.1 (for stan- 
dard needlets also for B = [1.2, 1.3, 1.4, 1.45, 1.5, 1.55, 1.7, 1.8, 1.9]), as well as the Bernstein 
needlet parameter k G [1,2,4,5] and the Mexican needlet parameter p G [1,2,3]. It will be 
shown that the critical angle is inversely proportional to i*. 

We therefore make a fit to the critical angle (which is measured from where the mask 
ends) of the form: 

ecru = f3/i\ (10) 



^This is more relevant for Mexican needlets, since they are not compactly supported as standard ones, as 
can be seen in the two top figures in figures |3l 
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where /3 is a parameter to be evaluated. Note that /3 depends on the threshold r and 
B, i.e. /3 = (]{B,t) (for Mexican needlets we found that /3 = P{t,p) with no B dependence 
but with dependence on p). We found that for standarcj^ needlets /3 can be fitted to obtain 
the critical angle for arbitrary thresholds r and values of B. The fit is of the form: 

m-) = C^(5)-((^^^ + C(r)), (11) 

where the 3 parameters a, ^ and ( depend only on the threshold r and the type of mask, 
while the parameter C depends only on B. Their functional form will be presented in the 
respective appendixes for the galactic and point source masks. For the mexican needlets, 
since they are largely independent of B, but depend on p, it is possible to make a fit of the 
fornj^: 

/3(r,p) = Ci(p)-r^^(p). (12) 
We will now discuss the details for each mask. 



4.1. Galactic cut 

In order not to get influenced by mask asymmetries, we use a symmetric mask extending 
15° on each side of the equator. Note that in section 1473) we will show that this modelling of 
a galactic cut works also for more realistic irregular galactic cuts like in the kq85-mask used 
by the WMAP-team or a Planck-like galactic cut. 

In appendix |D] we show the form of the hyperbolic fits we have obtained for the critical 
angles with this cut and compare the fits with actual calculated angles. Here we will only 
show the results for a selected number of cases. For standard needlets we will use 4 different 
values of B spanning from 1.1 to 2; for Mexican needlets we find no dependence of B and 
will show results for p = 1,2, 3. As the Bernstein needlets appear to be similar to the 
standard needlets, we do not show any plots for these, but show the fits for some choices of 



^The same holds for the Bernstein needlets (except for the additional slight dependence on the parameter 
k), but since their behaviour is very close to the one of the standard needlet, we refrained to make runs for 
many values of B in order to be able to make reasonable fits. Instead we report the values of /3 in table ([3]). 

■^Except for the two point source masks with p = 1, where a fit of another form was necessary. See 
appendix |E1 
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parameters in the appendix. The hyperbohc fits are not vahd at the lowest multipoles; the 
exact multipoles for which the fits are valid will be discussed in the appendix. 

In figure HI we show in the upper row the critical angles for standard needlets. In the 
left plot which shows the lowest multipoles, we see that for the most stringent threshold 
T = 10"'^ where basically no influence of the mask is accepted, the lowest scales i* < 40 no 
part of the sky can be accepted for any value of B. If a 10% influence is accepted, then 
an extension of up to 20 degrees may be sufficient for lower multipoles. We can clearly see 
how the localization is improved with higher values of B. The right plot shows the higher 
multipoles. For the most stringent threshold, an extension of 3-5 degrees is still necessary 
at i = 1000 whereas for a 10% influence, less than a degree is sufficient. 

In the lower part of the figure, we see the corresponding plots for the Mexican needlets. 
As expected, the higher the p, the worse the localization properties, but in all cases the 
Mexican needlets outperform the standard needlets in localization. Even for the most strin- 
gent threshold, for multipoles i < 20 there are still parts of the sky which may be used. At 
i = 1000, extensions less than one degree are accepted for all thresholds. 

4.2. Influence of the point source hole 

To mask point sources, the WMAP team uses circular holes of radius 0.6° (WMAP- 
hole); this radius corresponds to 2.5 times the beam FWHM (Full Width at Half Maximum) 
of 14'. The Planck channels with highest resolution will have beams of FWHM 5' and we will 
thus simulate Planck point sources holes with a radius of 0.21° (Planck-hole) corresponding 
to 2.5 times 5'. To test the influence of such holes on the needlet coefficients, we have placed 
a hole at the north pole. Note that in section H^Sl we will show that this modelling of a single 
point source mask works well also for more realistic point source masks, including many 
point-source holes. 

In appendix [E] we show the form of the hyperbolic fits we have obtained for the critical 
angles with point source holes and compare the fits with actual calculated angles. In the 
appendix, we also present the multipole ranges over which the fits are applicable. Again, we 
will here only show results for a range of cases. 

In figure Owe show the critical angles for the WMAP- hole and in for the Planck- hole. 
We see that for standard needlets, for i* < 300, the critical angle is normally at least a few 
degrees for all thresholds. For a mask with many point source holes, this practically means 
that for these scales, the whole map is affected. It is thus only for the smaller scales that 
a mask extension has any meaning. For Mexican needlets, this limit is for lower multipoles 
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Fig. 4. — Plots of some critical angles for the standard needlets (upper plots) and Mexican 
needlets (lower plots) for the galactic cut. The critical angle is defined to be at the border 
of the mask. A range of values in B, p and the threshold r has been chosen. The plots on 
the left show the smaller multipoles whereas the plots on the right are zoomed in on the 
larger multipoles. 

i ~ 100 — 200, depending on the threshold. For Mexican p = 1 and threshold of 10%, the 
mask extension is still less than one degree. 

For the largest scales, the hole does not significantly influence the needlet coefficients; 
the influence is so small that the above model breaks down for small multipoles. In table 
[T]we show the £*-ranges which are unaffected by the hole for different acceptance levels of 
influence. This is important information when working with the largest scales: when the 
multipoles are below these limits, the point source holes may be ignored. These limits are 
based on a threshold Tum defined as 

Tum= rrf^|l-C~(^)|, 
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Fig. 5. — Plots of some critical angles for the standard needlets (upper plots) and Mexican 
needlets (lower plots) for the WMAP point source hole with radius 0.6°. The critical angle 
is defined to be at the border of the mask. A range of values in B, p and the threshold r 
has been chosen. The plots on the left show the smaller multipoles whereas the plots on the 
right are zoomed in on the larger multipoles. 

which gives the total fraction of influence outside the hole. All needlet scales corresponding 
to £*'s below the limit given in the table have a total fraction of influence smaller than the 
given Tiim- For instance if a 1% (integrated) influence from the hole is accepted, then all 
scales corresponding to i* less than the one given in the table for = 0.01 may be used. 



4.3. Critical angle for more realistic masks 



As pointed out in the 2 preceding subsections the masks we used for finding the critical 
angles 6crit and fits to them were idealized. In reality one is interested in more complicated 
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Fit of critical angie 6, PLANCK-hoie, standard needlets Fit of critical angle 8, PLANCK-hole, standard needlets, high /* 




10 100 100 1000 

/• /■ 

Fig. 6. — Plots of some critical angles for the standard needlets (upper plots) and Mexican 
needlets (lower plots) for the Planck point source hole with radius 0.21°. The critical angle 
is defined to be at the border of the mask. A range of values in B, p and the threshold r 
has been chosen. The plots on the left show the smaller multipoles whereas the plots on the 
right are zoomed in on the larger multipoles. 

types of masks. In this section we show that the results obtained are also valid with good 
precision for more general masks. 

To investigate this issue we implemented the following procedure: 

• Choose a needlet, a frequency j, a threshold r, a value of B and obtain the correspond- 
ing critical angle 9crit using the fits given in the plots above or in the equations in the 
appendices. 

• Choose a realistic mask M and extend it with 6crit along the border to obtain Mmod- 

• Run simulations to obtain the correlation (7™""™ for the realistic mask and define an 
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extended realistic mask Mreai by setting to zero all points of the map ""^ where 
the value of the correlation is less than defined by the given thresholqj. 

• If the 6crit obtained from the idealistic mask is the correct angle also for a more realistic 
mask, the two extended masks Mmod and Mreai should be similar. In order to measure 
how well the angle works for the realistic mask, we take the difference Mmod ~ Mj-eai- 

We did this for the following masks: 

1. WMAP point source mask 

2. WMAP KQ85 galactic cut 

3. For a Planck-like galactic cut, we used the small KP12 mask used for the WMAP 1 
year release. 

We have tested several values of r, B, j and p for both galactic and point source masks 
and find excellent agreement between the extended masks; this suggests that the idealized 
models above work very well for realistic circumstances. For the lowest multipoles, the 
difference grows. Only at the lowest values of j where the mask extensions anyway are so 
large that they are of little use for realistic analysis, does the model start to be significantly 
different from the realistic mask extension. We will show three examples here at three 
different scales and for three different masks. 

In figure [7] we show an example of the difference Mmod — Mreai for WMAP-point sources 
for standard needlets with high multipole. The grey area shows the original mask and the 
blue/red points show where they differ at the border of the mask extension. One can see 
that the idealized WMAP-hole works well in zones with many holes. As expected, the model 
works less well when the distance between the point sources are comparable to the critical 
angle Ocrit, but if 6crit is used as a mask extension, this is of little significance as the full 
space between the sources will be masked anyway. 

In figures |8] we show selected results for the the galactic masks. Again we see that the 
difference between the critical angle and the realistic angle is always smaller than the mask 
extension, but that the difference starts to increase for the very smallest multipoles. 



''More precisely, for a given threshold r, set to zero all points of the map C™ where C™ "™ < 
CJ^~""^{6cr) with (7™-"™ being the correlation of the idealistic masks from the preceding two subsections. 
This takes into account that we define the critical angle via the total influence (oc |1 — C'™~"™(6')|(i0) 

and not only via CJ'-''"\ecrit) > 1 - t. 
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Fig. 7. — Check of generality for WMAP point sources holes for standard needlets with 
B=1.6 at multipole I* = 768 and with threshold r = 0.01. The plot shows the difference 
between M^od ~ ^reai- The light grey area shows the original mask. 

5. Harmonic localization properties 

In order to study the harmonic localization properties we looked at o"^ defined the 
following way: 



where is the needlet function in harmonic space. In the calculations below we have used 
'^max — 4500 (and Ngide — 4096). We remind that £* was formed as a weighted average of i 
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WMAP kq85 galactic cut, standard, B=1.6, I* =300 




1.0 



PLANCK-IIke galactic cut, standard, B=1.6, I* =46 




1.0 



Fig. 8. — Check of generality for KQ85 galactic cut (above) and "Planck" galactic cut for 
standard needlets with B=1.6 the former at multipole £* — 300 and with threshold r = 0.01, 
the latter at multipole i* — 46 and with threshold r = 0.1. The plot shows the difference 
between Mj^g^ — M^eai- The hght grey area shows the original mask. 

with hi as the weight. The expression for cr^ is formed as the weighted average of {I — £*Y 
thus being a measure of the width of h^. Since the Mexican needlets do not have a finite 
support, we expect that their harmonic localization is worse than the one from the other 
two needlets; this is indeed shown to be the case in the following. 
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We find that one can write 
meaning that 

0crit<^h = constant. 

We can thus, as expected, immediately deduce that the better the pixel space localization, the 
worse the harmonic space localization. It also implies that when the pixel space localization 
is strongly or weakly dependent on some parameter B, p or k, this will also be the case for 
the harmonic space localization. 

We note that the localization of the Bernstein needlets in harmonic space for A; > 1 is 
better than the one for the standard needlets. Further, since ah of the Mexican needlets 
does not depend on B, but in the case of the Bernstein and standard needlets it does, when 
increasing B, it happens that the Mexican needlets for high values of p have better har- 
monic localization properties than the other needlets (according to this measure of harmonic 
localization). 



6. Correlation properties 

Adding on to the previous discussion, we will also calculate directly the correlation 
between needlet-coefiicients at a distance 6 on the sphere. Using the definitions of the 
needlet coefficients, it can easily be shown that for an isotropic field, the correlation function 
Cj{9), defined as the correlation between 13 jk and f3jk' when the angular distance between k 
and k' is 6, can be written as 



4n 

where Pe{cos6) is a Legendre polynomial, and Ci is the power spectrum of the harmonic 
coefficients. We have checked that this expression holds using simulations. In figure |9] we 
see a plot of the correlation function for some selected values of 6 for standard and Mexican 
needlets. Looking at similar plots for Bernstein needlets with k G {1,3,5} we found that 
there are only negligible differences in correlation properties between them and the standard 
needlets. 
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Fig. 9. — Correlation function for different values of 6. The solid line is standard needlets, 
the others are Mexican needlets with p = 1 (short dashes), p = 2 (dash dot) and p = 3 (dash 
dot dot dot). Upper left: 6 = 0.5°, upper right: 6 = 4°, lower left: 6 = 8° and lower right: 
e = 30°. 

7. Summary 

We have performed the following set of tests of localization properties of needlets when 
used with CMB data: 

• The distance of influence from a galactic cut on needlet coefficients was found from 
simulations and a range from conservative (allowing no influence) to less conservative 
'distances of safety' (expressed by the critical angle Ocru) from the mask. We express 
the resulting distances (critical angles Ocru) as functional fits of the form and have 
obtained analytical fits for /3 for standard and Mexican needlets (see expressions [TT] 
and [T2|) . Results for some values are shown in figure H] (standard and Mexican) and 
listed in tables [2] (Mexican) and |3] (Bernstein). 
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• In a similar manner, the 'distance of safety' from a point source hole (one WMAP- 
like, one Planck-like) was found; analytical fits for the resulting critical angles (see 
again expressions [TU|, [TT] and [T^D were obtained. Some values are shown in figures O 
(WMAP-hole) and [H] (Planck- hole) and listed in table (Mexican) and E] (Bernstein). 

• We have then also shown that our idealized masks and the results obtained from there 
can be used and are valid for more general and realistic masks. In particular, the 
results are found valid both for the WMAP KQ85 mask and for a smaller Planck-like 
galactic mask. The same holds for point source holes: The results obtained for one 
hole is shown to be correct also for point source masks with many holes. 

• We have also studied the correlation distance of needlet coefficients on the sphere. We 
use an analytical formula which we have found to agree well with simulations (equation 
[H]) to find the degree of correlation as a function of distance between coefficients. 
Results for some distances are shown in figure O 

• We have defined a measure cr^ (see equation [T^ of localization in harmonic space. 
This characterizes the number of multipoles contained in one single needlet scale j and 
is found to be proportional to i*. As a result (remembering that Ocrit is proportional 
to and that Ocrit depends heavily on real space localization), we find that, as 
expected, Ocrit x cxpix is roughly a constant such that when the localization in real space 
is improved, the localization in harmonic space worsens and vice versa. 

For a certain type of needlet, one can thus obtain the critical angle by (1) calculating 
i* for the given scale j, using equation [HI (2a) reading the critical angle off figures HI Oand 
E] OR (2b) use the analytic fits presented in appendix [D] and [El 

Our results for the different types of needlets can be summarized as follows: 

• Standard needlets: The standard needlets are much better localized in harmonic 
space than the Mexican needlets (but similar to Bernstein needlets). In fact, the 
contribution to a certain scale j is coming from a limited number of multipoles with no 
infiuence from multipoles outside this range. However, the penalty for high localization 
in harmonic space is that the real space localization is lower than for other needlets. 
The parameter B controls the localization properties: The higher the B, the higher 
is the real space localization and the larger are the multipole ranges included in each 
scale j. 

• Mexican needlets: The Mexican needlets are much better localized in real space 
than any of the other types. In harmonic space however, a large (in principle infinite) 
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multipole range contributes to each scale j. The Mexican needlets depend on the 
parameter p, larger p have worse real space localization properties but better harmonic 
space localization. For p = 3 — 4, the real space localization properties of the Mexican 
needlets approach those of standard needlets for high values of B. For p = 1 and high 
j, the Mexican needlets are almost identical to the Spherical Mexican Hat Wavelets. 

• Bernstein needlets: Bernstein needlets are in between the Mexican and standard 
needlets. Their real and harmonic space localization is similar to the standard needlets, 
sometimes it is slightly better, sometimes it is slightly worse. The Bernstein needlets 
depend both on the parameter B and (weakly) on a parameter k. As for the standard 
needlets, a higher value of B increases the number of multipoles included in each scale 
j and improves the real space localization. The parameter k may slightly improve 
or worsen the localization properties depending on the exact measure used. As an 
example we find that for the critical angle for a galactic cut, for k > 2 the angle for 
the less stringent thresholds increases for increasing k, whereas the angle for the more 
stringent thresholds decreases for increasing k. 



8. Acknowledgements 



We would lik e to thank S . K. N aess for useful discussions. We acknowledge the use 



of the HEALPix (iGorski et al.l |2005| ) package. FKH is thankful for an OYI grant from 
the Research Council of Norway. This research was partially supported by the ASI/INAF 
agreement 1/072/09/0 for the Planck LFI activity of Phase E2. We acknowledge the use of 
the NOTUR super computing facilities. We acknowledge the use of the Legacy Archive for 
Microwave Background Data Analysis (LAMBDA). Support for LAMBDA is provided by 
the NASA Office of Space Science. 



A. A numerical recipe for standard needlets 



We recall here b r iefly th e recipe for the standard needlet construct i on wh ich is advo- 
cated in lBaldi et a"l~l J2009ah : IPietrobon et al. I tood ): iMarinucci et al. I feoosh 



STEP 1: Construct the function 



fit) 



exp( 



1 



-j^j, -l<t<l 
0, otherwise 
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It is immediate to check that the function f{t) is C°° (i.e. smooth) and compactly supported 
in the interval (—1,1); 

STEP 2: Construct the function 

J\f{t)dt 



The function (f){u) is again C°°; it is moreover non-decreasing and normalized so that 0(— 1) = 
O,0(l) = l; 

STEP 3: Construct the function 

[ 1 if 0< t <^ 

^it)^{ 0(l-J^(t-i)) i< t <1 
[ if t >1 

Here wc arc simply implementing a change of variable so that the resulting function (f)(t) is 
constant on {0,B^^) and monotonically decreasing to zero in the interval {B^^, 1). Indeed 
it can be checked that 



25 .(^ 1 , _ f 1 for t = l 



B 



and 



B -V B' \-l for t = l 



^(;^) = 0(1) = 1, 



STEP 4: Construct 



(^(1) = 0(-1)=O; 



and for h{^) we take the positive root. In view of all the above four steps, we see that 
h{^) e C°° with this construction method. 



B. A numerical recipe for Bernstein Needlets 

The steps for constructing a compactly supported function h{^) having a bounded 
derivatives (i.e. 6({) e C", a < k + 1/2) can be described as follows: 

STEP 1: Define 

B-l/B 
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where B > 1. 

STEP 2: Define polynomials (of Bernstein) 

Example: for n — 1 
ioT n — 2 

B^'\t) = (1 - tf , B['\t) = 2t{l - t) , Bi'\t) = , 

for n = 3 

B^'\t) = (1 - t)' , Bi'\t) = 3^(1 - tr , Bi'\t) = 3t\l - 1) , (t) = . 
STEP 3: Define polynomials 

i=0 

where 

Q = 1 , for i = 1, k, and Cj = otherwise. 

Example: for A; = 1 

1 

P3{t) = Y.Bf\t) = {l-tf + 3t{l-t)\ 

i=0 

2 

p,it) = Yl = (1 - + - 1)' + io^'(i - ■ 

Note that 

P3(0)=P5(0) = 1, 

and 

P3(l) = P3(0)=0, 

p',{t) = 30t\l-t)^ 
Therefore the derivatives of p^ have value zero up to order 2 at points and 1. In general 

P2ti(l)=P2ti(0)=0forr = l,...,fc. 
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STEP 4: Define the function 



if ^ > 5 
STEP 5: Finally, we have 

, otherwise 

In view of the examples given in this construction, we can also easily see that the 
constant Cm is in the order of M*^ (B — l/B)*^, which is rather small compared with the 
one could obtain from the recipe of standard needlets. 



C. A comparison between SMHW and Mexican Needlets for p = 1 



As mentioned earlier, it is suggested from results in ( IGeller &: Mayeli I ( l2009al )) that 
Mexican needlets in the special case where p = 1 provide asymptotically a very good ap- 
proximation to the widely popular Spherical Mexican Hat Wavelets (SMHW), which have 
been used in many papers on CMB analysis. More precisely, the discretized form of the 
SMHW can be written as 



'2' 



[2 - |,]e-^/--^ 



where the coordinates y = 2 tan | follows from the stereographic projection on the tangent 
plane in each point of the sphere; here we take 6 = 6jk [x) := d{x,^jk)- Now write 



by following the arguments in lGeller fc Mayeli I (l2009al ) and developing their bounds further, 
it can be argued that 



|*,-,(^; B-^) - K,k^P,k;i {0)\ = B~^0 (min {O^B^^l}) 



(CI) 



for some suitable normalization constant Kjk > 0. Equation fICip suggests that the numerical 
results in this paper in the special case where p = 1 can be used as a guidance for the 
asymptotic theory of random SMHW coefficients at high frequencies j. 



- 30 - 



D. The hyperbolic fits for the galactic cut 

In this section we will list the hyperbolic fits to the critical angle for the galactic cut. 
We use fits of the form 



where the 4 parameters of the fit C, a, ^ and ( will be listed in this section. For the Mexican 
we use fits of the form: 

/3(r,p)=Ci(p)-r^^(^). (D2) 

For standard needlets we obtained fits of /3 for the values B = [1.1,1.2,1.3,1.4,1.45, 
1.5,1.55,1.6,1.7,1.8,1.9,2] as weU as the thresholds r e [0.1, 0.01, 10"^, 10-^ 10"^]. We then 
fitted these values of /3 to the functional form in eqEJand obtained the following fits for a, 
^ and C, 

• C{B) = -0.3375 + 4.98 

• ^('^) = 0-89 independent of threshold; 

• a(r) = 0.059 ■ r-^-^^; 

• C(r) = 0.44 ■ r-O-32. 



We have tested several values of B and r which were not used for obtaining the above 
fits, and also in these cases the critical angles were found to agree very well with the critical 
angles obtained with the above fits. 

In order to show how well the fits work, in figure [10] we show some typical curves for 
the critical angle as a function of i* for different thresholds. We show both hyperbolic fits 
(of the form 6' = /?/£*) as well as the actually calculated points for the critical angles. 

For the Mexican needlets we obtain the following values of the constants Ci^2 for eqJT^ 



• ^1(1) = 2.2, Ci{2) = 2.4 and Ci{3) = 2.7 

• ^2(1) = -0.10, C2(2) = -0.14 and ^2(3) = -0.15 



We have thus obtained analytical formulae for the critical angle for standard and Mexi- 
can needlets. The Bernstein needlets have a very similar behaviour to the standard needlets 
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and we do not make a full fitting procedure for these, but results for some values of B and 
T are listed in table [31 

In tables [3] and H] we show the lower limit of i* for which the hyperbolic fit for the 
critical angle is applicable. This lower limit is due to the fact that at lower £*'s the galactic 
cut has influence over the whole sphere. In order to be more precise, we consider the mask 
to have influence over the whole sphere if the critical angle is bigger than TO'ifl. One would 
expect that this lower limit is higher for the more conservative thresholds, which is indeed 
the case. 



Fit of critical angle 9, galactic cut 




10 100 1000 



Fig. 10. — Plot of some critical angles of the standard and Mexican needlets for the galactic 
cut for different values of B, p and thresholds. The lines are the hyperbolic fits and points 
are measured values. 



""The maximal angle possible with our definition of the critical angle being zero at the border of the mask 
is 75°. 
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E. Fits for point source holes 

The functional form of the parameters entering equation [TT] to fit /3(i?,r) are: 

• Both holes: 

- C{B) = -0.3375 + 4.98 

• WMAP-hole 

- ^(t) = 0.89 independent of threshold; 

_ ^(^\ — ^^\fir] „n+V, f(^\ — -0.77 ^^^^ ( 1 Gogio(T)+2.43)^ Y 

a[r) - 10 ^ ^ with /(r) - -^_=exp j, 

- C(r) = 0.29 • 

• PLANCK-hole 

- Ht) = 0.57 independent of threshold; 

- a{r) = 10/M with fir) = ^^^exp (i ('°^io(y)+4.36)^ >| . 

- C(r) = 0.18 • 

which we obtained by fitting to the same set of values for B and r G [10~^, 3 -10"^, 10~^, 10~^] 
as for the galactic cut. 

For the Mexican needlets we obtain the following fit: 

• WMAP-hole 

- p=l : /3(r) = 4.3 ■ lO-^^^^ with 

fir) = -0.19 ■ (logio(r))' - 1.0 ■ logio(r) - 1.3; 

- p = 2 : /3(r) = 1.7 ■ r^O-^O; 

- p = 3 : Pir) = 2.1 ■ r'^-^^ 

• PLANCK-hole 

- p=l : /3(r) = 4.3 ■ lO^W with 

fir) = -0.16 ■ (logio(r))' - 0.89 ■ logio(r) - 1.2; 
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- p = 2 : /3(r) = 1.5 ■ r"' 

- p = 3: I3{t) = 1.9 ■ 



As for the galactic cut, we did not make a fit of the parameter /3 for the Bernstein needlets, 
as they are very similar to the standard ones. Instead, to give at least a glimpse, we report 
the values for some values of B in table E] only for the WMAP-hole. In table Owe show some 
values of /3 for the Mexican needlets (which can also be obtained via the formula). 

In figure [TT] we show some examples of the fits and actually calculated points. We show 
some of the worst as well as some of the better fits. In tables O [6] and [7] we show the lower 
limits on multipoles for which the models work. 



Fit of critical angle 9, WMAP and PLANCK hole 




100 1000 

z* 



Fig. 11. — Plot of some critical angles of the standard and Mexican needlets for WMAP-hole 
and Planck- hole for different values of B, p and thresholds. The lines are the hyperbolic fit, 
the points are measured values. 
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Table 1: Holes: Multipole ranges unaffected by the holes for different acceptance levels of 
influence, where 'W means WMAP-hole and 'P' Planck-hole. 55/60* means that the range 
is 55 for all k except k — 5. Hyphen means that no simulations were run for those values of 
the parameters. 
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Table 2: Galactic cut: Values of /3 for Mexican needlets for different thresholds. Below the 
lower limit on £* the mask has influence over the whole sphere. NB: /3 is given in radians. 
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Table 3: Galactic cut: Values of (3 for Bernstein needlets for different thresholds. The lower 
limit on I* is due to the fact that below those £*'s the mask has influence over the whole 
sphere. We list also the values of (3 of the standard needlets at the same 5's in order to 
show the small difference. NB: (3 is given in radians. 
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Table 4: Galactic cut: Lower limit on the (* for standard needlets. This limit on f is due to 
the fact that below those £*'s the mask has influence over the whole sphere. Hyphen means 
that there were not enough points which did not feel the influence over the whole sphere in 
order to obtain a trustable fit. 
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Table 5: Holes: Values of /3 for Mexican needlets for the different thresholds and the 2 holes. 
The lower limit on £* corresponds to the limit from where our model works. NB: /3 is given 
in radians. 
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Table 6: WMAP-hole: Values of ^ for standard and Bernstein needlets for different thresh- 
olds. The lower limit on i* corresponds to the limit from where our model works. We list 
also the values of /3 of the standard needlets at the same 5's in order to show the small 
difference. NB: /3 is given in radians. 
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Table 7: Holes: Lower limit on the ^* for standard needlets, where 'W stands for WMAP- 
hole and 'P' for the Planck-hole. This limit on ^* corresponds to the limit from where our 
model works. Hyphen means that we did not run simulations for that choice of parameters. 



